Transverse laser cooling of a thermal atomic beam of dysprosium 
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A thermal atomic beam of dysprosium (Dy) atoms is cooled using the 4/ 10 6s 2 (J = 8) — s> 
4/ 10 6s6p (J = 9) transition at 421 nm. The cooling is done via a standing light wave orthogo- 
nal to the atomic beam. Efficient transverse cooling to the Doppler limit is demonstrated for all 
observable isotopes of dysprosium. Branching ratios to metastable states are demonstrated to be 
< 5 x 10~ 4 . A scheme for enhancement of the nonzero-nuclear-spin-isotope cooling, as well as a 
method for direct identification of possible trap states, is proposed. 

PACS numbers: 06.20.Jr, 32.30.Jc, 37.10.De, 37.10.Vz 



I. INTRODUCTION 



Laser cooling and trapping of atoms has been a pro- 
lific area of research, with applications to the study of 
novel quantum phases 0, [J, matter-wave inter ferom- 
etry [If, quantum information and precision mea- 
surements and tests of fundamental physics Atoms 
amenable to direct laser cooling have recently expanded 
outside of the alkalis and noble gases to include elements 
such as chromium, indium, ytterbium and mercury ■ 
The lanthanidcs are of particular interest. The generally 
large ground-state magnetic moments of these elements 
present opportunities for new studies in areas such as de- 
terministic single atom sources (Toj . degenerate dipolar 
Fermi gases and quantum computing (l2| . While 
the complex electronic structures and numerous possi- 
ble decay channels for excited states do not make most 
lanthanides obvious candidates for laser cooling, some 
progress has been made. In erbium, the large magnetic 
moment allows for magnetic confinement in the trap re- 
gion while excited states recycle to the ground state, per- 
mitting laser cooling to low temperatures fl3f . A success- 
ful effort to form a Dy magneto-optical trap (MOT) has 
demonstrated that the same principle works for dyspro- 
sium [bij ]. 

In this work, we explore transverse laser cooling of dys- 
prosium (Dy, Z=66) in an atomic beam. In addition to 
sharing the largest ground-state magnetic moment in the 
periodic table with terbium (/i ~ 10 hb), Dy contains a 
pair of nearly degenerate, op pos ite-parity states the have 
been investigated in detail [15(- These levels have re- 
cently been used in a search for a variation of the fine- 



structure constant that does not rely on optical-frequency 
combs [lj| [l7| ■ The implementation of laser cooling in 
this experiment as a control over atomic-beam parame- 
ters is a crucial aspect of systematics management fl8| . 
Previously our group reported on the manipulation of a 
thermal beam of Dy atoms using a bright, visible transi- 
tion in Dy at 421 nm [r| (which was the first report of 
laser cooling in Dy) and the spectroscopy of this transi- 
tion (20j . In this work, we demonstrate successful trans- 
verse laser cooling of a thermal beam for all isotopes of 
dysprosium without repumping, propose a method for 
the simultaneous cooling and optical pumping within the 
hyperfinc manifold of the odd-neutron-number isotopes 
with non-zero nuclear spin, and show that the maxi- 
mum possible branching ratio of the cooling transition 
to metastable states is less than 5 x 10~ 4 . 



II. EXPERIMENT 

A. Dysprosium 

Dysprosium has seven stable isotopes, ranging from 
A = 156 — > 164, including two isotopes, 161 ' 163 Dy, with 
nuclear spin / = 5/2. Cooling is accomplished using the 
4/ 10 6s 2 (J = 8) — > 4/ 10 6s6p (J = 9) transition at 421 
nm [20l ]. Cooling of the odd-neutron- number isotopes 
was accomplished using the only closed transition in the 
hyperfine manifold from F = 10.5 -> F' = 11.5. The 



TABLE I: Parameters of the cooling transition 
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Parameter 






Value 


Wavelength 


A c 




421.291 nm 


Lifetime Exc. State [21j 


T 




4.8 ns 


Recoil Velocity 


V r — 


hk/M 


0.6 cm/s 


Doppler Limit 


Vd = 


y/h/2Mr 


19.8 cm/s 


Saturation Intensity 


1 sat 


= hcn/2r\ 3 c 


60 mW/cm 2 



2 



relevant cooling parameters of this system are given in 
Table H 



B. Apparatus 

This work was carried out using a beam source pre- 
viously used in a search for a variation of the fine- 
structure constant, a search for parity nonconservation 
(PNC), and rf-spectroscopy of nearly degenerate levels 
in Dy H3, ill HI- A detailed description of the 
atomic-beam source is given in Ref. [23j | . The beam is 
produced by an effusive oven with a multislit nozzle ar- 
ray operating at ~ 1500 K. The oven consists of a molyb- 
denum tube containing dysprosium metal, and is sur- 
rounded by resistive heaters made from tantalum wire 
enclosed inside alumina ceramic tubes. The atoms in the 
beam have a mean longitudinal velocity of ~ 3xl0 4 cm/s. 
The ID transverse velocity profile was measured via flu- 
orescence spectroscopy to be Gaussian (Eqn. Q}, with 
(T v = 1.3 x 10 3 cm/s. 

To generate the 421-nm light, 12 W of Ar-ion laser light 
(Coherent Innova 400) was used to pump a TkSapphirc 
ring laser (Coherent 899), producing up to 650 mW of 
842- nm light. After diagnostics and passing through an 
optical fiber, approximately 300 mW of 842-nm light en- 
tered a resonant bow-tie cavity containing a 1x2x10 
mm 3 periodically poled potassium titanyl phosphate 
(PPKTP) crystal (Raicol Crystals Ltd.) located at the 
22-/im beam waist of the cavity. The output power of the 
frequency doubler was consistently 50-60 mW. Outputs 
powers of >100 mW were observed, but at powers above 
60 mW, gray tracking in the crystal caused rapid degra- 
dation of the conversion efficiency [24| . It was determined 
experimentally that laser-cooling did not noticeably im- 
prove above ~ 30 mW of light power, so the cavity was 
operated at sufficiently lower power to eliminate notice- 
able degradation. The 421-nm light passed through a 
quarter- and a half-wave plate to adjust the polariza- 
tion to vertical and a 4x cylindrical telescope to create a 
0.5 mm x 1.5 mm rectangular beam profile . The 421- 
nm beam entered the atomic-beam apparatus and inter- 
sected the atoms at a ninety-degree angle in two regions 
(see Fig. [T}. The total length of the cooling region was 
approximately 3 cm, giving a total interaction time for 
an atom in the beam on the order of 100 /is. Assuming 
a cycle time of 4r ~ 20 ns for a saturation parameter of 
~ 1 [IH, this allows for 5 x 10 3 cycles. 

To reduce long term drifts of the laser frequency, 
the 421-nm light was scanned over the desired reso- 
nance while monitoring fluorescence with a photomul- 
tiplier tube (PMT) oriented at 90 degrees with respect 
to the laser and atomic beams. Once at the desired fre- 
quency, the TkSapph laser was stabilized using a scan- 
ning transfer-cavity lock [2f| • In this control system, light 
at 842 nm from the TkSapph laser is injected into a con- 
tinuously scanning Fabry-Pcrot (FP) cavity along with 
light from a mctrological HcNc (Spectra-Physics 117A) 
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FIG. 1: (Color online) Simplified layout of the experimen- 
tal system. Light from two Ar-ion lasers pumped a Coherent 
899 Ti:Sapphire and Coherent 699 dye laser respectively. The 
output of the Ti:Sapph laser was split, with 10% of power 
going to diagnostics and the frequency stabilization system. 
The rest was fiber coupled to the frequency doubling system, 
the output of which was expanded and polarized before en- 
tering the vacuum chamber. Light from the dye laser was 
fiber coupled to the beam apparatus, where it passed through 
polarization optics and an iris diaphragm before entering the 
chamber. 



with a frequency drift rate of < 0.5 MHz/hr. The trans- 
mitted light is detected with a photodiodc and sent to a 
peak-detection and timing circuit that measures the time 
delay between HeNe and Ti:Sapph transmission peaks. 
These time data are sent to an IsoPod microcontroller 
(New Micros, Inc.) which then adjusts the frequency 
of the TkSapph to keep the time delay constant. The 
frequency drift of the TkSapph stabilized this way was 
measured to be < 2 MHz/hr fl]|. 

After passing through the cooling region, the atoms 
traversed ~ 15 cm to the back of the apparatus for de- 
tection. The Doppler profile of the atomic beam was 
characterized by scanning a probe-laser frequency over 
the 4/ 10 6s 2 (J = 8) -> 4/ 9 5d6s 2 (J = 7) transition at 
658 nm. For a Gaussian velocity distribution described 
by 



P(v) = 



1 



2tt<t v 



(1) 



with characteristic half- width 1/e 2 velocity 2a v , the pure 
Doppler spectral profile is also Gaussian with the rela- 
tionship <7f = ka v /2n, where k = 2tt/X. The Doppler 
cooling limit, Vd = 19.8 cm/s (see Tabic IJ), corresponds 
to a spectral profile with a width cr Vn = 300 kHz at 658 
nm. Thus, the 658-nm transition is an ideal probe of 
the atoms' Doppler profile, as the natural width of the 
transition, 120 kHz [21j, is considerably smaller than the 
Doppler limited width. The sensitivity in our experiment 
was ultimately limited by the probe-laser linewidth; ef- 
fects such as magnetic-field broadening and power broad- 
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FIG. 2: (Color online) Example probe spectrum of the 164 Dy 
658-nm transition with 40 mW of cooling light. The un- 
cooled Doppler profile is shown by the unfilled diamonds. 
Optical pumping in the cooling region to magnetic sublevels 
with larger Clebsch-Gordan coefficients accounts for the non- 
conservation of area under the curve. The solid red line is a 
nonlinear least-squares fit to the data. The inset shows the 
narrow feature with the broad background subtracted and a 
single Voigt profile fit to the remaining peak. 



ening were controlled and maintained at a smaller level. 

The 658-nm light was produced using a Coherent 699 
dye laser pumped by 5.5 W of Ar-ion laser light. Af- 
ter passing through an optical fiber and diagnostics, 
~ 50 fiW of vertically polarized light in a 4-mm diameter 
beam entered the apparatus and intersected the atomic 
beam at a ninety-degree angle. Fluorescence from the 
probe transition was detected using a PMT oriented per- 
pendicular to the laser and atomic beams, with a 660- 
nm central wavelength interference filter (10-nm band- 
width) on the detection port. The dye-laser frequency 
was scanned continuously over the transition, while the 
fluorescence signal from the PMT and the signal from a 
reference photodiode (PD) were simultaneously acquired 
using an Agilent DSO5014A oscilloscope. Normalizing 
the fluorescence signal with the reference PD signal en- 
sured that fluctuations in the probe laser power did not 
contribute to features in the Doppler profile. 



C. Results 

A scan over the 658-nm 164 Dy resonance is shown in 
Fig. [21 The broad uncooled Doppler profile has a Gaus- 
sian profile with characteristic width, Of i = 20 MHz, 
corresponding to a Vi = 1.3 x 10 3 cm/s. At low tempera- 
tures the measured lineshape is expected to be described 
by a Voigt profile, a convolution of the Doppler spectrum 
of the atoms and the Lorentzian lineshape of the probe 
laser, given by 
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FIG. 3: (Color online) Probe spectra of the a) 163 Dy and b) 
161 Dy 658-nm, F = 10.5 — > 9.5 transitions with ~ 40 mW 
of cooling light. The unfilled diamonds show the uncooled 
Doppler profile, and the solid red line is a nonlinear least- 
squares fit with a Voigt profile used for the narrow feature. 
The slope on the right side of a) is from the edge of the 162 Dy 
transition. 



V(J) = (2) 
r 1 ( (f-f') 2 \ 1 7l/2 , 

where cr/ c and 7l are the Gaussian width of the cooled 
thermal distribution and the full-width half-maximum 
(FWHM) of the Lorentzian describing the laser spec- 
trum, respectively. For the best cooling, a good fit was 
achieved using a combination of a Gaussian and a Voigt 
profile: a Gaussian profile for the uncooled background 
and a Voigt profile for the narrow feature (see inset in 
Fig. [2j). The uncooled background is due to atoms outside 
the capture range of the cooling laser. The values for ct/ c 
and 7l obtained from a least-squares fit are 0.8(5) MHz 
and 4.2(7) MHz, respectively. The value for cr/ c gives 
<j Vc = 53(33) cm/s, which is consistent with the expected 
Doppler cooling limit of 19.8 cm/s (Table [I]). The value 
for 7^ is also consistent with a separate measurement of 
the laser linewidth (Sec. Ill D)) . 

Laser cooling of 161 Dy and 163 Dy, both of which have 
hyperfine structure, was accomplished by tuning the cool- 
ing laser onto the closed F = 10.5 —> 11.5 transition, 
where F is the total angular momentum (Fig. [3]). The 
Gaussian widths of the Voigt profiles fit to the narrow 
features of 163 Dy and 161 Dy spectra are 0.7(7) MHz and 
0.9(4) MHz, respectively, and agree with the Doppler 
cooling limit. 
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FIG. 4: (Color online) Power spectral density of frequency 
fluctuations of the dye laser used to observe the ID tempera- 
ture of the atomic beam. The solid red line is a least-squares 
fit to a function of the form given by Eqn. [3] 



The number of cycles can be estimated from the initial 
velocity, (a Vi = 1300 cm/s), the final velocity, and the 
recoil velocity given in Table U as N ~ (a Vi — a VJ )/v r « 
a Vi /v r ~ 2000. This allows us to constrain the branching 
ratio as 1 /N < 5 x 10~ 4 . This result is in agreement with 
a more sensitive measurement of branching ratios that 
has been performed in a MOT, yielding 7 x 10~ 6 [l4j |. 

D. Broadening Mechanisms 

Our measurement of the laser cooled transverse tem- 
perature of the atomic beam was ultimately limited by 
the linewidth of the dye laser. The linewidth was mea- 
sured by coupling the 658-nm light into a FP cavity sta- 
bilized to a metrological HeNe (Spectra-Physics 117A). 
The dye laser was tuned to the side of a transmission 
peak, where the slope is linear and the fluctuations in 
the transmitted power were recorded with a spectrum 
analyzer. The known free-spectral range (FSR) and mea- 
sured finesse of the cavity allowed us to to convert the 
voltage power spectral density (PSD) of the fluctuations 
to the frequency PSD [S v (f)]. The result of this mea- 
surement is shown in Fig. |4] A function of the form, 

SAf) = ^j L + S (3) 

was fit to the data, where 5*_i and So are the 1//- and 
white-noise contributions, respectively [23] ■ The least- 
squares fit values for S-i and So are 2.2(1) x 10 12 Hz and 
9.4(4) x 10 5 Hz respectively. Interpretation of the shape 
of the laser linewidth from this data is non-trivial, but 
the FWHM can be estimated by solving the equation [28| 

lrad 2 « ^ ^fidf, (4) 



where 2f^dB is the FWHM of the laser lineshape. Solv- 
ing Eqn. H for f 3dB we find a FWHM of 3.2(5) MHz, in 
agreement with the value obtained from a Voigt profile 
fit (Sec.lHCl). 

The most important remaining broadening mecha- 
nisms of the probe transition linewidth are magnetic-field 
broadening and power broadening. The magnetic-field 
broadening is given by 

Av B = \ {m Ju g u - m^gi^ixsB, (5) 

where g u = 1.26 and gi = 1.24 are the g- factors of the 
upper and lower energy levels and mj u and mj l are the 
magnetic sublevels involved. The worst-case scenario of 
driving m — > m ± 1 transitions gives a maximum broad- 
ening of ~ 1.9 MHz/G. This effect was minimized by 
coils installed around the probe region, which reduce the 
background magnetic field to ~ 1 mG. 

The other broadening mechanism of concern is power 
broadening. The saturation intensity, Is, of the probe 
transition is approximately 50 /iW/cm 2 , and the corre- 
sponding power broadened linewidth can be calculated 
using 



7p = 7o\/l + k, (6) 

where 70 is the natural linewidth and k = I /Is is the 
saturation parameter. With an approximate beam ra- 
dius of 2 mm, an input power of 50 /iW gives n ~ 8 
and 7p ~ 360 kHz, well below the linewidth of our laser. 
In addition, data were taken as a function of probe laser 
power to confirm the absence of measurable power broad- 
ening. 

E. Hyperflne Structure 

The dense structure of the 421-nm cooling transition 
gives rise to interesting optical pumping dynamics (spec- 
troscopic information on the 421-nm spectrum can be 
found in Ref. [H). In 161 Dy, all of the AF = transi- 
tions are unresolved due to the broad natural linewidth 
and Dopplcr width. Because the cooling transition 
is a A J = +1 transition in the uncoupled angular- 
momentum basis (J is the total electronic angular mo- 
mentum), driving all of the AF = transitions makes 
the ground state with the smallest value of F a dark 
state. In Fig. [5j the 421-nm laser has been tuned onto 
the cluster of AF = transitions while the probe laser 
is swept over the entire 658-nm spectrum, demonstrating 
optical pumping of the atoms into the F = 5.5 ground 
state. 

This suggests a convenient method for simultaneous 
cooling and optical pumping of 161 Dy. The 421-nm car- 
rier frequency can be tuned to drive the F = 5.5 — > 6.5 
transition, while a single modulation at ~ 1 GHz can be 
used to generate a sideband in the vicinity of the AF = 
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FIG. 5: (Color online) Example of the optical pumping dy- 
namics observed, a) A typical spectrum of the 658-nm tran- 
sition. For the odd-neutron number isotopes, the strongest 
AF = — 1 transitions are labeled with the ground state an- 
gular momentum, b) The same spectrum with the 421-nm 
light tuned into the vicinity of the AF = transitions. See 
Section Til El for details. 
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FIG. 6: (Color online) Diagram of the 163 Dy hyperfme struc- 
ture. The relative transition strengths ||dF|| 2 /||rfj|| 2 are 
shown alongside transition arrows, where \\d\\ is the reduced 
dipole matrix element normalized so that the F = 10.5 — » 
11.5 strength is equal to 1. The frequencies of the F — s- F + 1 
transitions relative to the F = 10.5 — > 11.5 transition were 
calculated from known hyperfme coefficients [2(| HiJ. The 
AF = — AJ transitions are negligibly weak and are not 
shown. 



transitions. The sideband will optically pump atoms to- 
wards the F = 5.5 ground state, while the carrier cools 
these atoms on the F = 5.5 — s- 6.5 transition. This has 
the potential to eliminate the multiple sidebands that 
would normally be required to cool the entire hyperfme 
manifold, as described in the next section. 



III. FUTURE WORK 

The extension of laser cooling to include all observable 
isotopes of Dy has been demonstrated, with cooling of the 
non-zero-nuclear-spin isotopes, 163 Dy and 161 Dy, accom- 
plished using the closed F = 10.5 — s- 11.5 transition. To 
improve the efficiency of cooling the odd-neutron number 
isotopes, we propose a scheme for cooling on the closed 
F = 10.5 — > 11.5 transition while simultaneously trans- 
ferring the population to the F = 10.5 state from some 
or all of the other hyperfme ground states. A diagram of 
the hyperfme manifold for 163 Dy is shown in Fig. [5J 

This can be accomplished by passing the 421-nm 
light throu gh a running-wave (non-resonant) electrooptic 
modulator [301 ] . The desired modulation frequencies can 
all be applied to the same crystal to generate sidebands 
at the appropriate frequencies. Atoms being driven on 
the F — > F + 1 transitions will eventually be pumped to 
the F = 10.5 ground state where they will continue to 
cycle. The effectiveness of this method will be limited 
by the power available to drive the modulator, however 



methods such as those described in Sec. Ill El may allow 
us to circumvent these limitations. 

Given the complex electronic structure of the dyspro- 
sium system, the existence of a closed optical transition 
is unlikely. A theoretical calculation was performed that 
estimated the total branching ratio to the most likely trap 
states at less than the 5 x 10~ 5 level [3l[ and a measure- 
ment in a MOT obtained a value of 7x 10" 6 [3 . Detailed 
knowledge of the dominant trap states will be important 
for optimization of MOTs and future work with laser 
cooled and trapped Dy. 

Our current set-up is insensitive to trap states for 
which the time of return to the ground state is less than 
the relatively long time-of-fiight, T to f ~ 500 /us, between 
the cooling and probe regions. We have developed a 
procedure to directly identify trap states. The 421-nm 
light will be tuned on resonance with the 164 Dy transi- 
tion while light from a pulsed dye laser operating between 
600 and 700 nm will be superimposed over the fluores- 
cence region. A PMT with a UV short-pass filter on the 
viewport will be oriented perpendicular to the atomic 
beam and the cooling laser while the dye laser is scanned 
from 600 to 700 nm. An increase in the fluorescence sig- 
nal from the PMT will be taken to indicate a transition 
from the trap state to a highly excited odd-parity state, 
the fluorescence from the ground state decay channel of 
this state being less than 300 nm. A diagram of likely 
repump wavelengths is presented in Fig. [7] 
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FIG. 7: Partial energy level diagram showing possible repump 
wavelengths from the likely trap states. The circled state is 
the upper level of the cooling transition. The trap states 
and corresponding repump wavelengths were identified in the 
NIST Atomic Spectroscopy database. 



cooling, but no narrowing of the line was observed due 
to the large probe laser linewidth. The remaining dom- 
inant broadening mechanisms of the probing method, 
background magnetic-field and power broadening, were 
estimated and controlled for. A plan for increasing the 
efficiency of cooling the odd-neutron number isotopes of 
Dy has been presented. The branching ratios of the cool- 
ing transition have been confirmed to be < 5 x 10~ 4 and 
a procedure for direct identification of the dominant trap 
states has been outlined. 

Future work includes using a femtosecond frequency 
comb as a probe of the Doppler profile of the atomic 
beam. The narrow width of the comb teeth should be 
sufficient to directly observe Doppler and sub-Dopplcr 
cooling. The ultimate goal of this work is the extension 
to 2D transverse cooling and longitudinal slowing in an 
experiment dedicated to a search for variation of the fine- 
structure constant. 



IV. DISCUSSION 

We have presented results on transverse cooling of a 
thermal dysprosium beam. A final laser linewidth lim- 
ited spectrum of a narrow probe transition in the cooled 
beam consistent with cooling to the Doppler limit was 
obtained. Cooling has been demonstrated for all de- 
tectable isotopes of Dy. A separate measurement of the 
laser linewidth was consistent with the hypothesis that 
the laser linewidth is the limiting factor in measuring the 
Doppler spectrum. An attempt was made at sub-Dopplcr 
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